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A single-party strategy in a multi-round quantum protocol can be implemented by sequential 
networks of quantum operations connected by internal memories. Here provide the most efficient 
realization in terms of computational-space resources. 

PACS numbers: 



Many results in Quantum Information and Quan- 
tum Estimation Theory 0,11] have been achieved through 
the general description of states and measurements in 
terms of density matrices and positive operator-valued 
measures (POVM's), respectively. The advantages of 
this formalism are evident in optimization tasks, like e.g. 
state discrimination, where one can look for the opti- 
mal measurement without considering the specific de- 
tails of the apparatus. Furthermore, the optimization of 
preparation/measurement devices is reduced to the op- 
timization of positive operators, for which many power- 
ful techniques are known. Similar advantages are pro- 
vided by the description of physical transformations as 
quantum channels (completely positive trace-preserving 
maps), which in turn can be represented by positive op- 
erators via the Choi-Jamiolkowski isomorphism f^]- 

The usage of the Ghoi-Jamiolkowski isomorphism is 
well established in quantum information theory Q 
since the early works on ancilla-assisted tomography 
0, @]. Recently, the Choi-Jamiolkowski representa- 
tion has been extended to more complex quantum de- 
vices, consisting of sequences of channels, quantum o 
erations and POVM's connected by internal wires 
11 1 . In particular, Ref. P considered the application 
of these sequential networks to represent single-party 
strategies in multi-round quantum games, while Refs. 
[lol 11 1 showed how these networks can implement a 
variety of higher-order quantum information processing 
tasks, such as transforming states into channels, chan- 
nels into channels, and even networks into networks. 
Refs. 
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also coined the name quantum combs for 
the Ghoi-Jamiolkowski operators associated to sequential 
networks, and developed a simple set of rules to describe 
the interlinking of networks in terms of the corresponding 
operators. In this framework, once a specific task is fixed 
(e.g. clonirig a channel [12] or estimating the POVM of 
a detector [13|) one can search for the quantum protocol 
that optimally realizes it. Having a simple description 
now becomes indispensable: since a quantum protocol is 
implemented by a complex network of devices, optimiz- 
ing each device separately is not a viable approach. In 
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FIG. 1: A quantum comb with slots. Information flows 
from left to right. The causal structure of the comb implies 
that the input system m cannot influence the output system 
n if m > n. 



the new framework, instead, the optimization of the pro- 
tocol is reduced to the optimization of a single positive 
operator subject to linear constraints. In the simplest 
cases the search can be also implemented automatically 
through matlab routines IJ, ll5| . 

Once the optimal Ghoi-Jamiolkowski operator has 
been found, however, one needs a way to unzip the in- 
formation contained in it and to find a physical imple- 
mentation of the network. In this Letter we solve this 
problem, presenting an automatic procedure that, given 
the Ghoi-Jamiolkowski operator of a quantum network, 
allows to construct a concrete implementation of it as 
a sequence of elementary devices. Among all possible 
implementations, the present one minimizes the compu- 
tational space, that is, at each step it uses the smallest 
possible dimension of the Hilbert spaces. Our procedure 
can be fully automatized in a computer software, accept- 
ing as an input the Ghoi-Jamiolkowski representation of 
the network and providing as an output the matrix rep- 
resentation of the operations that must be performed at 
each stage of the protocol. After the operations in the 
network have been determined one can look for a further 
decomposition of them into elementary gates, using e.g. 
the techniques of Refs. [Sill. 

We now review the basic concepts and results of the 
gen eral theory of quantum networks as presented in Refs. 
13, U I • The most general quantum device is a quantum 
circuit board, namely a network of quantum devices with 
open slots to which variable sub-circuits can be linked. 
By stretching and rearranging the internal wires of the 
network, we can give to each quantum circuit board the 
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FIG. 2: Linking of two combs. We identify the wires with the 
same label. 



shape of a comb, like in Fig. [T] The empty slots of the 
circuit board become the empty spaces between two teeth 
of the comb. Referring to Fig. [1] each wire is labeled with 
a natural number, which is even for the input wires and 
odd for the output ones; the corresponding Hilbert spaces 
are labelled in the same fashion (that is, the Hilbert space 
of the system represented by the wire i is denoted by 
Hi)). The ordering of the slots results from the causal 
ordering defined by the flow of quantum information from 
input to output; with our notation we have that input 
system in wire i can influence the output system in a 
wire j > i but not in a wire k < i. Two circuit boards Ci 
and 62 can be connected by linking some outputs of Ci 
with inputs of 62 , thus forming a new board 63 := 61*62. 
We adopt the convention that wires that are connected 
are identified by the same label (see Fig. [2]). 

In the following we will often use the isomorphism be- 
tween linear operators in Lin('H) and states inH ®7i: 
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where is a fixed orthonormal basis. 

The quantum comb C associated to a circuit board 6 
with N input/output systems is the Choi-Jamioikowski 
operator of the multipartite channel representing the 
input/output transformation that the board performs 
from states on 'Hin := ® j=o^ to states on 'Hout '■— 

^2i+i, T~in being the Hilbert space of the n-th sys- 
tem. A quantum comb is then a positive operator acting 
on 'Hout ® T-iin and it is defined as follows: 



Co 



:= (6®Ii„)(|/))((/|i„in) 



(1) 



(for clarity here we use the notation Hat = "Ho <8) V-b, 
Aab to mean A G hm{'Hab), to mean G Hb, and 
\A))ab to mean |^)) £ Hab)- It can be proved that the 
causal structure is equivalent to the recursive normaliza- 
tion condition 

Tr2fe_i[CW] =/2fc-2<»C('=~'^ fc = l,...,7V (2) 



where C(^) = C, C^o' = 1, CC^' e CiUo 



with 



^ink = (8>j=o "^2^ and Ho 



®j=o ^2j+i, is the comb 



of the reduced circuit C^'^^ obtained by discarding the last 
N -k teeth. 

The connection of two circuit boards is represented by 
the link product of the corresponding combs Ci and C2 , 
which is defined as Ci * C2 = Trjic [Cf C2] , Ojc denot- 
ing partial transposition over the Hilbert space JC of the 
connected systems (we identify with the same labels the 
Hilbert spaces of connected systems). 

One can wonder whether each positive operator which 
satisfies Eq. ([5]) corresponds to a sequential network of 
quantum channels. The answer is indeed positive, as 
shown in Refs. [§-ll| with the following Stinespring di- 
lation theorem: 

Theorem 1 Let C^^^ be a positive operator on 'Hout 
Hin, with Hin := <S>j'jo^ and Hout ■= ®^Jq'^2o+i- 
If C^'^^ satisfies Eq. then it is the Choi-Jamiolkowski 
operator of a sequential network given by the concatena- 
tion of N isometrics: for every state p G Lin('Hin) one 
has 

C W (p) - Tr^„ [y W • • • l/(i)py(i)t . . . y Wt] (3) 

where V'^'^^ is an isometry from 7^2^-2 T^Ak-i to 
^2fe-i (Xi Hyifc, and Ha^ an ancillary space, Haq — C 
(in Eq. (Q) we omitted the identity operators on the 
Hilbert spaces where the isometrics do not act). 

This result, however, provides little insight on how to 
construct the isometrics. We now give the explicit con- 
struction in terms of the Choi-Jamiolkowski operator in 
a way that can be automatically evaluated by a computer 
routine: 

Theorem 2 The minimal dimension of the ancilla space 
T-Lau *^ Theorem [7] is the dimension of the support of 
C^''\ A choice of isometrics V^''^ : H2k-2 ^Hau-i -> 
^2fe-i ®'HAk with minimal ancilla space is obtained by 
taking "HAk = Supp(C'-'^-'*), where * denotes the complex 
conjugation in the canonical basis, and 



l^))(2fe-l)(2fc-l)'T'(2fe-2)^(2fe-2)' 



(4) 



where T„_ 



Ei \i)m{i\n- 



has 



K and Eq. @ 

-1)-) = hk-2®lAk-i- 

Therefore, V'^^^ is an isometry. Now, define the 
isometry 

p^(fc) ^ which goes from 

T~Lini. to 'Houtfc ® 'H-Ak- By definition one has 



Proof. One 

(C('=-i)*)"^Tr2fe„i[CW*] 
yields yWt\/W = /2fc_2« 



Supp(C( 



(outfc) (outfc)' 71nfc^(infc)'J- 

However, according to Ref. [31, this is the minimal 
isometry of the channel C^'^^. Since the isometry is 
minimal, it is not possible to choose an ancillary space 
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FIG. 3: . The quantum circuit C, when Hnked with the uni- 
tary channel IX : Lin('Hi) Lin(^2), tries to reproduce the 
action of from input Ho to output Hs- 



smaller than 71 Ak- Finally, since C*^^^ is nothing but the 
channel associated to the network, Eq. ([3|) follows. ■ 

Theorem [2] implies Theorem [l] and provides a recipe 
for the concrete realization of the quantum network with 
minimal dimension of the ancillas at each step. The di- 
mension of the ancilla is the quantum "space" of the com- 
putational network. Note that sometimes the isometries 
V^''^ can act trivially on some subsystem, this resulting 
in further simplifications of the physical implementation. 

As an application of the methods outlined above we 
now consider the problem of finding the quantum net- 
work that realizes the optimal inversion of a unitary op- 
eration. Such a network consists of a circuit board 6 
with an empty slot to be linked to the unitary channel 
U{p) — U pU'^ sending states on Hi to states on H2. The 
resulting circuit C * U has to be as similar as possible to 
the channel 1X~^ (see Fig. ([3])). The quantum comb of 
e is C e LinCH32io), with ^ n2 ^ V-i no C, 
and, according to Eq. satisfies the normalization 



Tr3[q =/2®C(i), Tri[CW]=/o. 



(5) 



Choi operator of the unitary channel is \U)) {{U\2i and the 
link e*U gives the operator (([/* |2iC|[/*))2i G Lin(-H3o). 
To quantify the similarity between the channel 6 * U and 
the target we use the channel fidelity [l^: if A 
and B are two channels and A and B are their Choi- 
Jamiolkowski operators the channel fidelity T{A,B) is 
defined as f{d~^A,d~^B) where / is the state fidelity 
a) = I Tr ^/^fofhJaS^ . In our case we have 

F(C,Z^) - /(d-i(CHC/»((C/|2i),d-i|[/t))((C/t|3o) 



d2 



((C/^|30((C/*|2lC|C/t))30|C/*))2 



(6) 



We assume the unknown unitary \J randomly distributed 
according to the Haar measure of SU{d), and, as a figure 
of merit, we adopt the average of the gate fidelity: 

F = [ dUF{C,U) 

Jsu{d) 

1 



d2 



SU(d) 



dC/((C/^|30((t/*|2lC|C/t))30|C/*))21 (7) 



where dU denotes the invariant Haar measure. The fol- 
lowing lemma holds: 



Lemma 1 The operator C maximizing the fidelity ^ 
can be assumed without loss of generality to satisfy the 
commutation relation 

[c, U3(g>W2(g>Ui(g> Wo] = yv,w e SU{d) . (8) 

The proof consists in the standard averaging argu- 
ment: Let C be optimal. Then take its average C — 
JdUdW {U3(g)W2(g>Ui(8)Wo)CiU3(g>W2(8)Ui(g>WQ)^: 
it is immediate to see that C satisfies Eqs. ^ and (O, 
and has the same fidelity as C. 

Thanks to Schur's lemmas C can be decomposed as 



(9) 



where S = {+,—}, is the projector onto the symmet- 
ric/antisymmetric subspace of Hi €5 V-j , and a^'' ^ 
Vp, V. Moreover, using Eq. ^ the fidelity ([7]) becomes 



1 



^=^(ai30((/|2lC|/))30|/»21 

^^a'^^d,, d. -TriP"- 



d2 



(10) 



while the normalization ([5]) becomes X^^jes '^^"'^a' = 
l,Vi^ G S. The last equality implies the bound F — 
X^^eS ^'^'^d.fj. ^ 2/(^2^ which is achieved if and only if 
a^'^ — ^4^, that is, if and only if 
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20 



d^ (i_ 

dU \U^)){{U^30(^\U*)){{U*\21. 

SU{d) 



(11) 



We now use Theorem [5] to construct the optimal network 
from the quantum comb C. Since C(i) = d-iJio the first 
isometry is given by 

namely it consists in the preparation of the maximally 
entangled state while the input state is stored 

in a subsystem of the ancilla space Hai C Hyo' . 

The second isometry F^^) . -^^ (g) Hai ^ Tis (S) 'Ha2 is 
given by 



y(2) = (Vd/3^C*^)|/))33, ®T2. 



(12) 



Remarkably, this is the Stinespring isometry of a 
measure-and-prepare channel. Indeed, consider the chan- 
nel £{p) = Tr^jy^^VV^^^^I and the POVM 

Q. = {C*)-^U^)){{U^\3-o' ® \iJ)){{U\2'v{C*)-K (13) 
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FIG. 4: . Optimal circuit for the inversion of a unitary trans- 
formation, the input state in wire is stored in a quantum 
memory. The unitary U to be inverted is estimated and the 
inverted estimated unitary is applied to the input. 

which provides a resolution of the identity in 'Ha2 = 
Supp(C*) due to Eq. We then have 

= d j dUU\{U\2'vp\U))2'vU, (14) 

namely the channel £ can be implemented by measur- 
ing the POVM P^j = d\U)){{U\2'v on the Hilbert space 
7^2'!' and subsequently performing the unitary U'^ on 
Hoi. Therefore, the optimal network for the inversion of 
an unknown unitary channel corresponds to an "estimate 
and re-prepare" strategy: first the isometry V^^^ pro- 
vides the optimal input for the estimation of U (that is, 
the maximally entangled state d^^\I))iii), then, after the 
unknown unitary has been applied, the second channel £ 
performs the optimal POVM on the state d^^^Ujjiy and, 
depending on the estimate C/, applies the unitary on 
the input state stored in wire 0'. The physical implemen- 
tation involving measurement and classical feed-forward 
is an alternative to the coherent, fully quantum process- 
ing corresponding to the isometry T^^^-*. 

In conclusion, we provided a general method for re- 
covering all the isometries of a network from its Choi- 
Jamiolkowsky operator minimizing the comptational 
space. This result allows us to formulate an algorithm 
for designing optimal quantum networks for any desired 
task (e. g. cloning, discrimination, estimation): 

1. Choose a suitable figure of merit F for the task of 
interest. 

2. Find the positive operator C satisfying constraint 
in Eq. ^ and maximizing F . 

3. Set C(o) = 1 and Ia„ = 1. 

4. For k = 1 to k = N do the following 

(a) Calculate /^T ® CC') = Tr^[C], where % 
(Tr:j^) denotes the identity (partial trace) over 
all Hilbert spaces but % 



(b) Define V^*^) as in Theorem [2] 

5. The optimal network is given by the concatenation 
of the V^^'^^'s in Eq. © 

We applied the algorithm to design the optimal circuit 
for the inversion of a unitary transformation. It is worth 
noting that in general the numerical optimization of step 
2 can be challenging, and that it is typically convenient 
to exploit the symmetries of the problem to reduce the 
number of parameters, as we did here in our example. On 
the other hand, the remaining steps 3-5 — which represent 
the original result of the present Letter — can be easily 
programmed on a computer. 
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